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We determine the quantitative effect, in the maximum number of particles and other static ob- 
servables, due to small anharmonic terms added to the confining potential of an atomic condensed 
system with negative two-body interaction. As an example of how a cubic or quartic anharmonic 
term can affect the maximum number of particles, we consider the trap parameters and the results 
given by Roberts et al. [Phys. Rev. Lett. 86, 4211 (2001)]. However, this study can be easily 
transferred to other trap geometries to estimate anharmonic effects. 
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I. INTRODUCTION 

The experimental realization of Bosc-Einstein conden- 
sation (BECj in magnetically trapped weakly interact- 
ing atoms fl|-|| brought a lot of interest in theoretical 
studies of the properties of condensed atomic systems. 
For systems with negative two-body scattering length, re- 
cently it was reported in Ref. |)J a discrepancy between 
experimental results and theoretical predictions [||, in 
the maximum critical number of trapped atoms, N c . In 
Ref. Q|, it was shown that part of the discrepancy can 
be explained by taking into consideration the axially de- 
formed shape of the trap. The theoretical prediction of 
the critical number of atoms, for the cylindrical symme- 
try considered in Ref. Q, has to be adjusted to a number 
that is lower than the number given by the corresponding 
spherically symmetric trap. Still, this correction is not 
enough to obtain a result that is totally compatible with 
experimental values. So, it is relevant to look for other 
possible sources of the observed discrepancy, in this spe- 
cific for example, higher order non-linearities in 
the mean- field approximation J8| , or possible experimen- 
tal deviations not already taken into account. 

The parameter associated with the critical number N c , 
given in Ref. |J, is defined by 

k= , (!) 

yfi/ (muj) 

where uj = {oj x oj v oj z )^^ is the geometrical mean value of 
the trap frequencies, m the mass of the atomic species, 
and a the two-body scattering length of the particles 
in the condensate. With the assumption of a spherical 
symmetrical trap, k = 0.575 The critical number 
of atoms for Bose-Einstein condensates with a < have 
been investigated by the JILA group, taking into account 
experiments with 85 Rb [^), considering a wide tunning of 
the scattering length a, from positive to negative values, 
by means of Feshbach resonance j9). The experimental 
result for the number fc, k exp t = 0.459± 0.012 (statistical) 
±0.054 (systematic), corresponds to a deviation of about 
20% lower than the predicted spherical symmetrical re- 
sult. Part of this discrepancy was shown to be related to 



the cylindrical cigar-shaped symmetry employed in the 
JILA's experiment, where the frequencies corresponding 
to the directions z and r± = \J x 2 + y 2 are, respectively, 
given by lo z — 2n x 6.8 Hz and u>±_ — 2tt x 17.35 Hz. So, 
the parameter of anisotropy is given by A = lo z /loj_ = 
0.3919. The theoretical prediction, obtained in Ref. 0, 
with the correct cylindrical symmetry given in Ref. Q , is 
k = 0.55. In this case, the experimental result is smaller 
than the predicted value by about 16.5%. Such a result 
is still not enough to include the theoretical result within 
the error bars of the experimental result. 

As explained in Ref. [[l0|, in the JILA's experiments 
with 85 Rb, the trap is practically harmonic in the cen- 
tral region, for very low temperatures and for small size 
condensates. However, one should check how a devia- 
tion of the harmonic trap, outside of the central region, 
can affect the number of condensed atoms. In this case, 
we are considering the possible existence of imprecisions 
in the form of the confining potential, generated by the 
modified Ioffe-Pritchard design used (known as baseball 
trap). The modification of the ground state solution of 
the condensed state can alter, correspondingly, the ob- 
servables associated with it. 

In our present investigation we consider a deviation in 
the harmonic trap potential that is effective only outside 
the central part of the potential, with the addition of a 
term that is proportional to a cubic or quartic power of 
the distance. This work was first motivated by looking 
for a possible source of the observed discrepancy between 
theory and the JILA's experiment. However, we realize 
that, if exist any deviation from the harmonic trap in the 
experiments reported in Ref. ||, it should be of a very 
small factor considering the kind of trap design used by 
them pp[ . But, one should be aware also of other kind of 
trap arrangements in the experimental studies of BEC, 
where an investigation of possible effects in the observ- 
ablcs due to anharmonic terms in the interaction can be 
useful. In this perspective, our present study of the effect 
of anharmonic terms added to a harmonic trap interac- 
tion, is not restricted to the example that we are going to 
consider. We consider the trap parameters of the JILA's 
experiments || as an example, estimating deviations in 
N c due to anharmonicities, that can easily be extended 
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to other geometrical trap configurations, with the help 
of previous studies [Q. Presently, we are reporting nu- 
merical results obtained for N c when the trap deviates 
from the harmonic shape. We should also note that there 
exists a previous study considering the occurrence of an- 
harmonic terms in a time orbiting potential trap . 

As shown by our results, a steeper confinement will re- 
sult in a lower k, against the naive expectation. At first 
sight, a steeper confinement should increase the kinetic 
and potential energies, with stabilization for a larger (and 
not smaller) value of k. However, a steeper confinement 
also means that the relative distance between the atoms 
is less than in the case of the harmonic potential. There- 
fore, with the atoms feeling more attraction from the 
interaction term, this will result in collapse for smaller k. 

Next, we review the formalism used in the present ap- 
proach, followed by the main results and conclusions. 



and the trap potential, Ut ra p 
can be rewritten as 
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where the normalization of is given by 
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In the critical limit, we have n — n c 

The equation, corresponding to Eq. (^), for the chemi- 
cal potential f3 in dimensionless units (p = fnojS) is given 
by 
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II. MEAN FIELD APPROXIMATION 

The Gross-Pitaevskii equation (GPE) that describes 
the wave-function of the condensate \P in the mean-field 
approximation has the form 



ay 

lh -dt = 



2m 



Utra 



4nh 2 \a\ 



(2) 



where the potential Utrap = Ut rap (f;u! x ,co y ,co z ) is given 
by a modified harmonic oscillator trap, and the wave 
function ^> = $(r,t) is normalized to the number of 
atoms N. For the stationary solution of Eq. (|J), with 
\t(r, t) = exp(—ipt/h)y(r,0), where p is the chemical 
potential, we obtain 



*n- = p*. 



(3) 



Considering the symmetry used in the JILA's experi- 
ment, where the trap frequencies are given by 10 ±_ = lo x — 
<-Oy 7^ l> z , it is adequate to work in cylindrical coordinates: 



\J x 2 + y 2 and 6 = arctan ^— ^ 



(4) 



For the ground state of the condensate, we have \& = 
^(rj_, z,t). By using the trapping geometrical average 
1/3 

frequency to = (u> x ujyU} z ) and the oscillator length by 
l Q = y/h/(jnui), we obtain the following dimensionless 
coordinates and parameters: 

p=V2^, C=%/2^, T = ut; 



Jp = , UJ C = 

LU CO 



(•5) 



With the above dimensionless units, redefining 



As known jf2|, the equation (g) is valid in the mean- 
field approximation of the quantum many-body problem 
of a dilute gas, when the average inter particle distances 
are much larger than the absolute value of the scattering 
length and the wave-lengths are much larger than the 
average inter particle distance. 

The total energy of the system is given by 
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In order to analyze the effect of a deviation of the trap 
potential from the harmonic behavior, we consider two 
expressions for V(p, £), labeled by v = 1, 2 as follows: 

v {v) ( P , C) = \ Hp 2 (1 + 8 pP v ) + ^C 2 (1 + 5 C C)] ■ 

(ii) 

The distortion added to the harmonic potential is cubic 
when v — 1, and quartic when v = 2. In both cases, the 
magnitudes of the distortions in the directions p and £ 
are given by the parameters S p and 5^. 

With Eq. ([ll]) , in both cases one can observe that the 
interaction is approximately harmonic near the center of 
the trap (at p = and/or ( = 0). When p < 1, the 
quartic term keeps approximately the harmonic shape of 
the potential in a more effective way than the cubic term. 
In Fig. nL we can see how the shape of the harmonic trap 
(S p = 0) changes in the radial direction, at the position 
C = 0, when one adds a cubic (y = 1) or a quartic (v = 2) 
term in the potential. 
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FIG. 1. The anharmonic potentials, as given in Eq. (p"l|), 
with 8 P = 0.2 and £ = 0, for the cubic (u = 1, dashed-line) 
and quartic [v — 2, dotted-line) case, are compared with the 
harmonic potential (S p = 0, solid-line) All the units are di- 
mensionless, such that Utrap = (fiio)V and r± — p l /s/2. 
The values of the frequencies correspond to the ones given in 
Ref. [§]: lu z = 2?r x 6.8 Hz and ui ± = 2tt x 17.35 Hz. 



III. RESULTS OBTAINED 

In this section, we report our main results consider- 
ing the solutions for the chemical potential, obtained 
for Eq. (^|) with the anharmonic interaction given by 
Eq. (|ll|), for several combinations of the parameters S p 
and Sq, for v — 1 and v = 2. We have also considered the 



solutions for the total energy, which is given by Eq. (10). 
For the numerical solution of Eq. (Q), which is given 
in cylindrical symmetry, we have employed the Crank- 
Nicolson algorithm in two dimensions, employing the re- 
laxation method propagating in the imaginary time, as 
earlier described in |7|]. The spatial discretization used 
was 140 x 140, with cut-off p m ax — 7 and (max = 7. 
The time discretization was At = 0.005 and the relax- 
ation time used for obtaining the solutions of stationary 
ground state was r = 16. The results are more sensitive 
to the spatial grid spacing, but a lack of precision occurs 
in the third decimal figure of the values obtained in the 
numerical simulations. With these grids, we obtain pre- 
cise results up to the second decimal figure, confirmed by 
the convergence obtained with a refined spatial grid (160 
x 160). We obtain our results for the total energy and 
the chemical potential of the condensate, as functions of 
the number of particles, by considering the geometry of 
the trap given in JILA's experiment reported in Ref. || , 
where lo z = 27rx6.8 Hz and ui± = 2ttx 17.35 Hz. We con- 
sider a range of values for the parameters 5 P and 5^, that 
represents the corresponding anharmonicity. As shown, 
only positive anharmonic terms added to the original har- 
monic trap can make the theoretical prediction for the 
maximum critical number k become closer to the JILA's 
experimental result. This implies that one should check 
for possible deviation of the harmonic potential, outside 



the central region, that makes the trap more confining. 

In Table I, we present our results for the maximum 
critical number k, which is related to the critical number 
of atoms N c , as given by Eq. ([!]), considering several 
possible choices for the anharmonic parameters S p and 
<5<j, for both cases v = 1 (cubic) and v = 2 (quartic). 

In the Fig. 2, we show some results for the case 
with v = 2, where we have a quartic anharmonic term 
added to the original interaction, in both directions p 
and £, such that S p — S^. In this example, we are also 
considering the cigar-like symmetry of the trap used in 
Ref. [|[, withwx/Ws = 17.35/6.8. We present two frames: 
for the total energy E (lower frame), given in units of 
(NTilo) / (2n) and for the corresponding chemical poten- 
tial \i (upper frame), in units of Tilo, versus N\a\/l a . 

This picture does not change too much when we switch 
off the deviation in the (, direction, with <5^ = (the trap 
remains harmonic in the z— direction). For the case with 
v = 1, where we consider cubic anharmonic terms added 
to the original interaction, our results show that for both 
5 Z = S p and 5^=0 there are no significant qualitative 
differences in the observables plotted and the values are 
approximately equal to the quartic case. The behavior 
is similar for all the cases considered, with the collapse 
of the condensate occurring for smaller critical number 
as we add larger positive deviation in the harmonic trap. 
In this respect, the addition of a quartic term {y = 2) in 
the potential, is more effective in keeping the harmonic 
shape near the center of the trap. 

TABLE I. Values of the maximum critical number k, for 
several combinations of the parameters S p and 5q , related to 
the anharmonic factors introduced in the interaction, as given 
in Eq. (p"l|), for v — 1 (cubic) and v = 2 (quartic). The 
cylindrical symmetry is the same given in Ref. M. 
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FIG. 2. Total energies E (lower frame), in units of 
(Nhuj) I (2n) [i? to tai = E(Nhcu)/(2n)], and chemical poten- 
tials fi (upper frame), in units of hui [fi = /3(hu))], versus 
N\a\/l a , for different values of 8 P — S^, as given inside the up- 
per frame. The deviations from the harmonic trap are given 
by quartic terms in both directions p and (,. Here, we are also 
using the trap considered in Ref. JH), with the asymmetry 
given by uj±/uj z = 2.5515. The corresponding critical num- 
bers are at the end of each curve, with the position of the 
experimental value k expt pointed out with vertical arrows. A 
corresponds to the sum of the experimental errors. 
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FIG. 3. Variations of the critical number of atoms, k, in the 
condensate, as functions of 5 P . The corresponding deviations 
in the £ direction, 8^, are given inside the figure. The position 
of the experimental k and the corresponding positive error bar 
(A) are also indicated. 

All the above results can be globally analized by ex- 
amining Fig. 3, in which we have the behavior of the 
maximum critical number of atoms in the condensate 



state, parametrized by k, as a function of the magni- 
tude of the trapping anharmonic potential, parametrized 
by S p . We display the effect of both cubic and quartic 
anharmonicity in confinement, with or without includ- 
ing the deviation in the z-direction. As we can see, the 
higher is the anharmonic parameter, the smaller is the 
critical number of condensed atoms. As one can observe, 
in order to obtain the theoretical results for k inside the 
region covered by the experimental error bars, one needs 
deviations of about 10% to 20% from the harmonic trap, 
at distances of the order of the oscillator length. 



IV. CONCLUSIONS 

In summary, we have solved numerically the station- 
ary Gross-Pitaevskii equation in cylindrical symmetry, 
for condensed systems with attractive two-body interac- 
tion, with anharmonic trapping potentials, either with 
cubic and quartic deviation from the harmonic oscilla- 
tor. We did the present study in the perspective to ob- 
serve the real effect on N c of anharmonicities that can 
occur in the usual confining traps that have been exten- 
sively employed in the experiments with atomic BEC. We 
took the harmonic trap parameters considered in Ref. [^) , 
to exemplify how much should be the strength of cubic 
or quartic anharmonic terms in the trap interaction to 
present a sizable change in the critical parameter fc; in 
the perspective of a partial explanation of the observed 
experimental deviations from the previous theoretical es- 
timates. The present results show that only with a de- 
viation of about 10% from the harmonic potential, at 
distances of the order of the oscillator length, it becomes 
possible to explain theoretically the experimental results 
reported in within the error bars. However, this is 
quite large strength for the anharmonic term, that is not 
supported by the detailed analysis and description of the 
experiments given in Ref. |l(| . 

In a more realistic perspective, our present work is pre- 
dicting the quantitative effect in the critical number of 
atoms due to anharmonicities that can occur in a trap. 
For the non-harmonic terms we have considered cubic or 
quartic ones. In cases of other shapes for the anharmonic 
terms, one can also make qualitative predictions based on 
the present results. We should point out that our study 
is not restricted to the example (trap frequencies and 
symmetry) that we have used; it can be easily extended 
to other trap arrangements, and can be useful in experi- 
mental analysis of condensates with attractive two-body 
interactions. In this respect, we should note that it has 
been studied and reported the possible occurrence of an- 
harmonicities in time orbiting potential traps [Q| . Our 
work is considering realistic situations, that are close to 
the actual experiments. So, our aim in the present work, 
was to report the effect of a deviation in the harmonic 
trap that increases as we go to regions outside the center 
of the system. 
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